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1 . SCOPE 


document describes the conversion of the CLASSY clustering/ 
classification program from the UNIVAC EXEC 2 in JSC Bldg 12 
to the IBM 370/148 computer at U\RS, West Lafayette, Indiana. 

Th0t converted program is written exclusively in Fortran IV-G. 

The conversion approach has been similar to that followed in 
the EOD-LARSYS system conversion. 

Several program enhancements and/or changes have been incorporated 
as a result of information made available during the conversion 
process. 
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2 . APPLICABLE DOCUMENTS 


• TIRF 77-0072 

• "As-Built" Design Specification, "CLASSY Program Modification," 
TIRF 77-0055, JSC-13986, LEC-12185, April 1978. 

• "Program Documentation for Modification to the CLASSY Program," 
JSC-12602, LEC-10481, April 1977. 

• Technical Memorandum, "Pinal Acceptance Test Plan for the 

EOD-LARSYS conversion," Ref: 646C-13, March 1978. 


3. SYSTEM DESCRIPTION 


3-1 HARDWARE DESCRIPTION 

The CLASY clustering program, as nsoclifled, is operational on the Uni vac 1110 
(or 1108) under the EXECS operating system. The program utilizes the Univac 
Fortran V compiler (the original program utilizes the Univac reentrant Fortran 
compiler, RFOR), the Univac assembler (for assembly routines FREE, GET 
and LOCK), and the Univac system random file access routines RINIT, RREAD, 
and RWRITE. 

3.2 SOFTWARE DESCRIPTION 

The CLASY system of subprograms was originated by Dr, Michael Rassbach, a post- 
doctoral contractor for NASA-JSC, Earth Observations Division (TF). 

The purpose for his development of CLASY was the iriiplementation of an iterative 
statistical clustering algorithm which had theoretical promise for application 
to classification of earth resources (image) data acquired from the LANDSAT 
(formerly ERTS) satellite. 

These modifications to CLASY are (1) to improve the program execution time, 
especially input/output overhead, in order to make it feasible to test and 
evaluate the CLASY program, (2) to implement a randomized Input data 
scrambling technique which was obtained from Rice University by the Earth 
Observations Division's Research, Test and Evaluation (RT&E) Branch (TF3). 

The scrambling of inpub data vectors is a necessity in CLASY due to the 
algorithms; sensitivity to correlated data values and (3) to allow users to 
specify more than 4 channels for classification. 

The CLASY system of subprograms consists of the main driver program, CLASY, 
and 56 subprograms, not including the Univac system routines utilized by 
the program. Four of the subrprograms (LOCK, GET, FREE, and BYTRAN) were 
originally programmed Univac assembly language (SLEUTH II), the remaining 
subprograms and the main program, CLASY, are converted to the Univac Fortran 
V language from the former (original) Univac reentrant Fortran language, RFOR. 


The driver program for the clustering system is CLASY. The data handling 
subprograms for the system are READTP and STATIS which were reprogrammed to 
implement the required data input improvement and to implement the randomized 
data scrambling technique. READTP provides the data setup for acquisition 
on demand by the iterative statistical subprogranu STATIS. STATIS initiates 
the clustering procedure, operating on one pixel (data vector) at a time 
in setting up clusters and making the cluster split/combine decisions. Each 
pixel is examined 10 times by STATIS during the clustering procedure, 

To implement the required modifications and the modifications suggested by 
the program originator to improve the program's reliability, the following 
routines in the CLASY clustering system received changes: C1.ASY, CLASYl 
(deleted), ADJUST, LOCK (formerly LOG), CLPRM, (formerly CLPR) CLDUMP, MISH, 
SETUP9, CMBKIO, ELIM, MULTI, SEPER, TR, STATIS, CLUST. 

The modified subprograms are discussed in the order of their use by the CLASY 
system, 

The overall CLASY system is flowcharted in Appendix A. Listings of the 
modified routines are shown in Appendix B. Sample output from the CLASY 
system is shown in Appendix C. 

3.2.1 SOFTWARE COMPONENT NO. 1 (CLASY) 

3. 2. 1.1 Linkages 

CLASY is the driver program of the CLASY clustering system. CLASY calls 
SETUP9, READTP, MULTI, and CLUSMP. 

3.2. 1.2 Interfaces 

The common blocks INFORM, CLUSTR, CLUS, MISC, and STPAR and calling arguments 
are used in the program CLASY as interfaces with other routines in the 
clustering system. 
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3. 2. 1.3 Inputs 

CUSY calls SETUP9, which reads the Input supervisor (control) cards. The 
supervisor cards and their functions are described In the discussion of 
SETUP9 (section 3. 2. 2. 6). 

The required Input to the CLASY program consists of one tape (or file) con- 
taining the multichannel Image data and the special -format card input, 

The Image data tape (file) Is presumed to be In either of two specific formats-- 
elther "LARSYS II" format or "UNIVERSAL" format. The tape (file) reading 
program In CLASY* TAPERD, accepts either of these formats and sel f-determlnes 
the correct method of reading the data. 

3. 2. 1.4 Outputs 

The output by CLASY is all line-printer output. Interim printout of statistical 
parameters and diagnostic data Is provided during the Iterative cl us ter- forming 
process . 

The final output Is a "map", with a symbolic representation of area clustered, 
with each pixel of the area classified using the statistics (mean and covariance) 
from final cluster set determined by CLASY. The symbols on the "map" represent 
the cluster (-class) which Is the most likely parent distribution for the 
given pixel. The "map" Is output by subprogram CLUSMP. 

Sample output Is shown In Appendix C. 

3. 2. 1.5 Storage Requirement 

Storage used: Code = Data = 42112g 

3. 2. 1.6 Description 

CLASY Is the driver program for the clustering routines. It was rewritten to 

(1) enable CLASY1 to retrieve large blocks of data (dimensioned ARRAY 

(20 ood)) from drum^ (2) to pre-calculate the amount of data RREAD (drum read 
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subroutine) will read from drum, (3) to wake DATAB (the array containing 
the scrambled data) and PV (the array passed to the clustering routine, jTATIS) 
to reside in the same locations as the large data array, "ARRAY" and (4) to 
selectively skip calculations on clusters which have subclusters. All other 
logic in the program remains the same as in the pre“modified version of 
CLASY. 

3. 2. 1.7 Flowchart 
See Appendix A. 

3. 2. 1.8 Listings 

See Appendix B for program. 

3. 2. 1.9 Restrictions 

The known restrictions inherent in the program are (1) the program will not 
successfully execute with only one channel, (2) a data vector containing a 
zero value in the channel of interest will cause an error termination of the 
program's execution, (3) the size of the original image data set read from 
the input tape (or file) and placed on drum must be containable in 1,310,717 
locations of drum storage available to the random access routines (RINIT, 

RREAD, and RWRITE). 

3.2.2 SOFTWARE COMPONENT NO 2 (SETUP9) 

3. 2. 2.1 Li Pkage 

SETUP9 is called from CLASY. SETUP9 calls NXTCHR and NUMBER, which are entry 
points in subroutine FIND. 

3. 2. 2. 2 Interface 

Interface is accomplished through calling arguments and the following common 
blocks: INFORM, SUPCUM, and CLUSTR. 


3.2.2. 3 Inputs 
See Appendix C. 
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3. 2. 2. 7 Flowchart 


See Appendix A. 

3.2 .2 ,8 Listings 

•See Appendix B for the modified program listing, 

3.2.3 SOFTWARE COMPONENT NO. 3 (READTP) 

3. 2. 3.1 Linkage 

READTP is called from CLASY. READTP calls RREAD, RWRITE, CMERR, and ZOR 
(Function ZOR is the random number generator used in the data scrambling 
technique), RINIT, TAPHDR, LAREAD, FLDINT, LINERD, FDLINT and ERTRAN. 

3. 2. 3. 2 Interface 

Interface is accomplished through calling arguments and the following common 
blocks: INFORM, CLUSTR, CLUS, MISC, and STPAR. 

3.2. 3. 3 Inputs 

Image data tape described in 3, 2. 1.3 
3.2.3. Output 

READTP outputs the following error messages: 

3. 2. 3. 5 Storage Requirement 

Storage used: Code = 1652g Data = 427^ 

3.2. 3.6 Description 

READTP performs the input image data-handling function for the CLASY 
clustering system and makes the image data available on disk to the iterative 
statistical subprogram, STATIS. The original image data from the area on 
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the input tape (file) which has been designated by the input field-definition 
card(s) is preston?d on drum as one continuous block of data. The Univac 
ranv-n file access routines — RINIT, RWRITE — are utilized to place the 
input data on the drum as the data is read from the input tape (file) by the 
TAPERD subprogram. 

READTP precomputes the base addresses for three data-buffering arrays in 
core-storage. One half of data-buffering array* ARRAY is used to retrieve 
a block of original image data from the prestored drum. 

The other half of array, ARRAY is used to contain a set of integers — 1,2,3,..,,N 
where N»the number of data vectors in the original image array, Upon completion 
the first half of the array will contain a set of data vectors obtained from 
the original image set, but stored such that each data vector's original 
spatial location is randomly rearranged. The data scrambling technique 
utilized in rearranging the data vectors was obtained from Rice University 
via RT&E and is the required modification implemented by this modification to 
the CLASY clustering system. The implementation of the randomized data 
scrambling technique is performed as follows: 

(1) Given an array of original image data vectors, and an array, of integers -- 
l,2,3n...,N with N - the number of image data (2) Scramble the 

(2) Scramble the elements of A 

a. Obtain a random number, Z^. , from the uniform random number generator; 

Z^.=Z0R(0), 0.<Z^<1.0 

b. Multiply the random number by N, the largest integer in A; 

. IX. =N X Z. + 1 

c. Using IX^. as an index, scramble the integers in A as follows: 
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3. 2. 2. 4 Output 

SETUP9 prints out a summary of the input to CLASY and also prints an error 
message on the line printer, if an invalid input card is detected If an 
error is detected, SBTUP9 prints the following message "INVALID INPUT CARD- 
IGNORED", processing continues, 

3 . 2 . 2 . 5 Storage Requi rement 

Storage used: Code « 2/^jj Data » 157^^ 

3. 2. 2. 6 Description 

SETUP9 reads and analyzes all cards input to the CLASY program. SETUP 9 was 
modified to add "0" (zero) as a new symbol to the symbol array "SMBLS". 

The following control cards are input to the modified CLASY program, to be 
analysed by SETUP9. In all cards, the "keyword" begins in card column 1, 
and any parameters on the card are placed from card columns 11 through 72, 
inclusive. 

1. "CHANNEL" CARD (i.e., "CHANNEL 1,5,9,13") 

The "CHANNEL" card specifies the channel numbers to be used in clustering the 
multi -channel data vectors. At present the maximum number of channels 
allowed to specify is eight. The identification "CHANNEL" starts in column 1, 
and the actual channel numbers, separated by commas, start in card column 
11, and must be terminated by column 72. 

2. "PRINT" CARO (i.e., "PRINT 1,3,3") 

The "PRINT" card specifies how to print the cluster map. The identification 
"PRINT" starts in column 1, the actual print parameters start in card column 11, 
separated by commas, and ending by column 72. 

3. "HEDl" card 
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4. "KED2" card 

These two cards may be used to specify any arbitrary heading for the printer 
output, including the cluster map. Any alphanumeric characters put into 
card columns 11 -72 of these two cards will be output as a page heading. 

5. "NPOS" card 

This card, used previously in the unmodified CLASY program, specified the 
number of positions to skip to read a pixel point (i.e. to "scramble" the 
data). This card is not used in the modified CLASY program. 

6. "NPTS" card 

This card, used previously in the unmodified CLASY program, specified the 
number of pixels to retrieve from the data set at each point. This card 
is not used in the modified CLASY program. (NPTS was used in the original 
programs data scrambling technique). 

7. "DATE" card 

This card is used to specify the date or any eight characters. Will be 
printed at the upper right hand corner of each page of printer output. 

8. "COMMENT" card 

The "COMMENT" card is equivalent in use and format with the "HEDl" and 
"HED2" cards, described above. 

9. "*END*" card 

This card specifies the end of all supervisor (control) card (described 
above) input to CLASY. This card is a mandatory input to CLASY, to 
initiate the clustering process. 








TEMP = 

= IX^ 

Aix. = TEMP 

d. Execute the «bove procedure (a-c) N times, with 

i = N, N-1, N-2, .... 1 

e. Create the new (scrambled) set of data vectors, A g, 
as follows; 

for i " 1 , 2 , 3 , . . . , N 
As^ = A(j), where j = Ai 


The scrambled data set, is made available in large blocks for retrieval on 
demand from the statistical clustering subprogram, STATIS. The scrambling 
of data for STATIS is necessitated by the clustering algorithm's sensitt ty 
to correlated data. 

The data buffering technique as described accomplished the first objective 
of the modifications to the CLASY clustering system— namely, the improvement 
of program execution time. 

3.2. 3.7 Flowchart 
See Appendix A. 

3.2. 3. 8 Listings 

See Appendix B for program listing. 

3.2.4 SOFTWARE COMPONENT NO. 4 (MULTI) 
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3.2.4. 1 Linkage 


MULTI 1s called from CLASY. MULTI calls DATFIX. ALFREE, CLINIT, STATIS 
and CLDUMP. 

3. 2. 4. 2 Interface 

Interface Is accomplished through calling arguments and the following 
common blocks J CLUS, MISC, STPAR, INFORM and CLUSTR. 

3. 2. 4. 3 Inputs 
None. 

3. 2. 4. 4 Output 
None. 


Storage Requirements 

Storage used: Code =^05 Data = H_ 

o 8 


3. 2. 4. 6 


Description 


MULTI calls the routines to initialize the clustering algorithm. MULTI 
was modified to make use of unused areas in core. ARRAY (EXARAY) was pre- 
viously dimensioned but never used. A method was devised to make use of 
this unused core such that the PV ana.^ and DATAB array (array utilized in 
STATIS and the array containing sorauMed data vectors) be made to utilize 
the same area of core (i.e., the PV and DATAB arrays were made equivalent 
storage areas, in ARRAY (EXARAY)). 

Flowchart 
See Appendix A. 

3'2.4.8 Listings 

See Appendix B for program listing. 

3.2.5 SOFTWARE' COMPONENT NO. 5 (STATIS) 

3. 2. 5.1 Linkages 

STATIS is called by MULTI. STATIS calls DISC, CLASY2, CORECT, DOTSQ, 

VPV, VMTV, MPVS, ADJUST, CLDUMP, and EXP. 

3-2-5.2 Interface 

Interface is accomplished through calling arguments and the following 
common blocks: CLUS, MISC, STPAR, CLUSTR, and RAND. 


S' 
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3, 2. 5. 3 Inputs 


None. 

3. 2. 5. 4 Outputs 

STATIS outputs two warning messages. They are: "***WARNIN6 ON THE 
INDEX(KL)=", "**SUSPECTED BAD DATA POINT -STATIS**IDO'=_.ROOT„, 
VECTOR " 


3. 2. 5. 5 Storage Requirements 

Storage used: Code = 1305g Data = 306^ 

3. 2. 5. 6 Description 

STATIS takes each input data vector and classifies it on a fractional, 
probabilistic basis. It then updates the various statistical parameters 
associated with the classes (clusters) indicated and checks to see if any 
of these classes is potentially two. Those which are will be referred 
to the' routine "SPLIT". The one modification made to STATIS was the re- 
dimensioning of the PV array, which contains the data vectors to be clustered. 

3.2. 5. 7 Flowchart 
See Appendix A. 

3. 2. 5. 8 Listings 

See Appendix B for program listing. 


S' 
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3,2.6 SOFTWARE COMPONENT NO,, 6 (ADJUST) 


3. 2. 6.1 Linkage 

Adjust 1s called from STATIS. ADJUST calls GET, TR, DOTSQ, SQMTX, MINV, 
UNIF, CLPR, TRIMTX, DENCAL, SPLIT, FREE, CLDUMP, SEPER, SUBLIM, ELIM, 
CORECT, JOIN, APRIOR, SQRT, ALDG, EXP, and XPRI. 

3-2. 6. 2 Interface 

Interface is accomplished through calling arguments and the following 
common blocks: CLUS, MISC, STPAR, CLUSTR, and JOINPR. 


3.2, 6. 3 Inputs 


None. 


3. 2. 6. 4 Outputs 


ADJUST prints out three brief messages concerning statistical Information 

and three error messages. They are; "ADJUST WEIGHT WAS SPFAC 

CHANGE "STATISTICS: TRACE__ SKEW__ KURT__ TESTS (SPLIT >0); 

"WHAVE SPLIT_ WEIGHT_ SUBS ", "W/OVOL ERROR IN ADJUST 

KL,W,NEW W,V0L ", "***EXTRAPOLATION PROBLEM IN ADJUST: ITER, 

INDEX(KL), VOLIN, OVOL, CV0L_„ ^", "LOG ERROR IN ADJUST: I, 

IM, KL, K/VRIN= " 

3. 2.6. 5 Storage Requirements 

Storage used: Code = “ ^^^8 
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3. 2. 6. 6 Description 

One modification made to ADJUST was to eliminate the use of subscripted sub- 
scripts for those arrays which previously presumed the RFOR (reentrant 
Fortran) compiler. The change was made to enable the routine to compile under 
the Un^vac 111-EXEC 8 Fortran V compiler. The second change made was that 
input to the "ALOG" routine is forced to be positive by use of the absolute 
value of the input being sent to "ALOG". The reason for the change was than 
an occasional negative value was being sent to "ALOG", causing an error 
termination. The routine previously referenced a dimensioned variable in 
several calling statements causing an error during program execution. These 
errors were corrected. Other modifications which were suggested by the 
CLASY program originator, Dr. Mike Rassbach, are changes in the calculation 
of EXF, WADJ, DCORR and ALINK. 

3. 2. 6. 7 Flowchart 
See Appendix A. 

3. 2. 6. 8 Listings 

See Appendix B for program listings. 

3.2,7 SOFTWARE COMPONENT NO. 7 (TR) 

3. 2. 7.1 Linkage 

TR is called from ADJUST. 

3. 2. 7. 2 Interface 

Interface is accomplished through calling arguments and the following common 
blocks: CLUS. MI SC, and STPAR. 

3. 2. 7. 3 Inputs 
None . 
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3.2. 7.4 Outputs 


None. 

3.2.7. 5 Storage Requirements 
Storage used; Code » 124g Data 35g 

3. 2. 7. 6 Descriptio n 

The modification made to TR was the elimination of subscripted subscripts, 
to enable this routine to be compiled under the Fortran V compiler on the 
Univac 1110 (EXEC 8 system). The subscripted subscript notation was a 
feature of the reentrant Fortran compiler RFOR, utilized by the program 
originator in development of CL.ASY. 

3. 2. 7. 7 Flowchart 
See Appendix A. 

3. 2. 7. 8 Listings 

See Appendix B for program listing. 

3.2.8 SOFTWARE COMPONENT NO. 8 (CLPRM) 

3.2.8. 1 Linkage 

CLPRM is called from CLDUMP, ADJUST, SEPER and JOIN. CLPRM calls GET, LOCK, 
SQMTX, MINV, and FREE. 

3. 2. 8. 2 Interface 

Interface is accomplished through calling arguments and the following common 
blocks: CLUS, MISC, and STPAR. 
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3.2.8. 3 inputs 


None. 


3. 2. 8. 4 Outputs 


CLPR prints out seven messages concerning statistical information. These 
messages are: "aUSTER__ INDEX_ PROPORTION^ SPLIT_ WEI6HT_ WAS__ 

A0JUST__ T0__ PROPORTION; PR0P__ CIN._ CTUT„ OLD PR0P_ CIN^,^ 0DEN_ 

DIFf'ER™, V0LUME__ R0UT_„ 0C0N__", "LOCATION_ LINK_ SUBS__ SUPER 

SYMBOL__", "NET PROB__ DIRECT^ CUMS "GUMS "MEAN 

"KURT(*W) ", "OLD COVARIANCE 

3. 2. 8. 5 Storage Requirements 


Storage used: Code *736*0, Data *364 g 


3. 2. 8.6 Description 

The modifications made to CLPR to form CLPRM were (1) the elimination of sub- 
scripted subscripts (2) a reference to one of CLASY's clustering routines, 
named "LOC", caused ambiguity, because there existed a Uni vac system routine 
with the same name. Therefore this reference and all other references to 
"LOC" were changed to "LOCK". The originator of CLASY, Dr. Michael Rassbach, 
provided changes to be made in CLPR. One modification suggested by 
Dr. Rassbach was the changing of a format statement so that the printer will 
skip a line before writing the statistical mean. 

3. 2. 8. 7 Flowchart 
See Appendix A. 

3. 2. 8. 8 Listings 

See Appendix B for program listing 
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3,3.9 SOFTWARE COMPONENT NO. 9 (CLDUMP) 


3.3.9. 1 Linkage 

CLDUMP is called from CLASYl, MULTI, CLASY2 and ADJUST. CLDUMP calls 
ISPLIT and CLPR. 

3. 3. 9. 2 Interface 

Interface is accomplished through calling arguments and the following common 
blocks: CLUS, MISC, STPAR, and CLUSTR. 

3.2.9. 3 Inputs 
None. 

3. 2. 9. 4 Outputs 

A header concerning the dump of the clusters and a debug statement giving 
the value of the number of the split cluster and the value of the print 
parameter. 

3. 2. 9. 5 Storage Requirements 

Storage used: Code = 153g, Data - 33g 
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3. 2. 9. 6 Description 


The modifications made to CLDUMP were the elimination of the logical "IF" 
statement that checked the value of ISPLIT and the print variable, PROUT, 
before the call is made to CLPR. A debug printout statement was added to 
allow printout of these variables for checkout purposes. 

3. 2. 9. 7 Flowchart 


See Appendix A. 

3. 2. 9. 8 Listings 

See Appendix B for program listing. 

3.2,10 SOFTWARE COMPONENT NO. 10 (ELIM) 

3.2.10.1 Linkage 

ELIM is called from ADJUST, ELIM calls SUBLIM and TRFREE, 

•S' 

3.2.10.2 Interface 

Interface is accomplished through calling arguments and the following 
common blocks: CLUS, MI SC, and STPAR. 

3.2.10.3 Inputs 
None. 
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3»2.10.4 Outputs 


EUM prints out two messages, one statement when a cluster has been eliminated 
and the other is an error message when a structural error has occurrjid. 

These messages a»'er ''#if#ELIMINATE_ LINK, LSUBS, LSUPER « j', 

"♦♦STRUCTURAL ERROR AT ELIM; KEL, KFAITH, KULD, INIT 

3.2.10.5 Storage Requirements 
Storage used: Code ® 176^, Data « 46g 

3.2.10.6 Description 

The modifications made to ELIM were the eliminations of subscripted sub- 
scripts in arrays to allow the routine to be compiled under the Fortran V 
compiler of the Univac 1110-EXEC 8 system. 

3.2.10.7 Flowchart 


See Appendix A. 

3,2.10.8 Listings 

See Appendix B for program listings. 

3.2.1! SOFTWARE COMPONENT NO. 11 (SEPER) 

3.2,11.1 Linkages 

SEPER is called from ADJUST. SEPER calls CLPR, DGNCAL, and FREE. 


S' 
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3.2.n,2 Interface 


Interface Is accomplished through the following common blocks: CLUS, MISC, 
and STPAR. 

3*2.11.3 Inputs 

None. 

3.2.11.4 O utpu ts 

SEPER prints one statement, which is printed whenever a cluster is split. 

This statement is : "##il»SEPARATE_ SUPER, SUBS SPFAC_" 

3.2.11.5 Storage Requirements 

Storage used: Code =237g Data ® 

3.2.11.6 Descriptions 

The modifications made to SEPER were the elimination of subscripted subscripts 
1n arrays, to enable the routine to be compiled under the Fortran V compiler 
of the Uni vac 1110-EXEC 8 system. 

3.2.11.7 Flowchart 
See Appendix A. 

3.2.11.8 Listings 

See Appendix B on program listing. 
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3.0.12 SOFTWARE COMPONENT NO. 12 (LOCK, fomerly "LOC) 


3.2.12.1 



LOCK is a function subprogram and is called from CLPR, and CLPRM, 

3.2.12.2 Interface 

Interface is accomplished through the function arguments. 

u 

3.2.12.3 Inputs 
None. 

3.2.12.4 Output 
None . 


3.2.12.5 Storage Requirements 
Storage used: Code IOiq 

3.2.12.6 Descriptions 

The modification made to LOCK (formerly "LOC") was the changing of the 
external reference name from 'LOC to 'LOCK*. The name 'LOC is also the 
name of a Fortran V system routine which caused ambiguity between the two 
routines during execution of CLASY. Changes were made to all other routines 
that referenced this routine, to change the reference to be "LOCK" rather 
than "LOC". 

3.2.12.7 Flowchart 
See Appendix A. 
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3.2,12.8 Listings 


See Appendix B for program listing. 

3.2.13 SOFTWARE COMPONENT NO. 13 (WISH)' 

3.2.13.1 Linkages 


The Fortran V procedure, MISH, is included in the following routines: 
ACOM, ADJUST, ALFREE, AMSQ, APF,I0R, CBLO, CLASY, CLDUMP, CLINIT, 

CLPRM, CLUSflP, CLUST, CORECT, DATFIX. DOTSQ, FlGROT, ELIM, ISPLIT, 

JOIN, MINV, MLT, MPVS, MTVEC, MULTI, MVEC, SEPER, SPLIT, SQMTX, STATIS, 
STOFLO, SUBLIM, TR, TRFREE, TRIMTX, VMTV, and VPV. 

3.2.13.2 Interface 


The Fortran V procedure, MISH, is used as the interface for the following 
common blocks: CLUS, MISC, and STPAR. 

3.2.13.3 Inputs 


None. 

3.2.13.4 Output 
None. 

3.2.13.5 Storage Requirements 

Storage used: None 
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3.2.13.6 Descriptions 

The modification made to the proc, MISH, was the changing of the symbol 
array 'NSYMB' dimension from 11 to "12" (i.e., NSYMB(12)). This was done 
because another symbol was added to the large symbol array, 'SYM*. 

3.2.13.7 Flowchart 
See Appendix A. 

3.2.13.8 Listings 

See Appendix B for program listing. 

3.2.14 SOFTWARE COMPONENT NO. 14 (CMBKIO) 

3.2.14.1 Linkages 

The Fortran V procedure, CMBKIO, is included in the following routines: 
ADJUST, CLASY, CLDUMP, CLINIT, CLUSMP, JOIN, MULTI, READTP, SETUP9, 
and STATIS. 

3.2.14.2 Interface 

The Fortran V procedure, CMBKIO, is used as the interface for the common 
block , CLUSTR. 

3.2.14.3 Inputs 
None. 

3.2.14.4 Output 


I 


t 






None. 


XX 


3.2.14.5 Storage Requirement 


Storage used: None 

3.2.14.6 Descri pti pn 

The modifications made to the Fort.r.m V procedure CMBKIO were as follows: (1) 

The changing of the parameter, MAXPQP, value from 60 to 61. This was done 
to increase the dimension of the symbol array to account for the addition of 
another symbol to the array. The symbol added to the array was the character 
'O' (zero). (2) The maximum number of channels was changed from 4 to 16. 

3.2.14.7 Flowchart 


See Appendix A. 

3.2.14.8 Listings 

See Appendix B for program listing. 

3.2.15 SOFTWARE COMPONENT NO. 15 (CLUST) 

3.2,15.1 Linkage 

CLUST is called from CLUSMP. CLUST calls ISPLIT, CORECT, DOTSQ, and EXP. 
3*2.15.2 Interface 

Interface is accomplished through calling arguments and the following common 
blocks: CLUS, MISC, STPAR, BIGCOM. 

3.2.15.3 Inputs 

None. 


3-24 


3.2.15.4 Outputs 


None. 

3.2.15.5 Storage Requirements 

Storage used; Code =332g , Data 8 

3.2.15.5 D pscrlptlo n 


3.2.15.7 Flowchart 
See Appendix A. 

3.2.15.8 Listing 

See Appendix B for program listing. 
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4.0 OPERATION 


The source code and relocatable elements for the CLASY system are on tape 
12667, located in Building 12, Johnson Space Center, Houston, Texas. The 
CLASY clustering program is operational on the Univac 1110 (or 1108) under 
the EXEC8 operating system. 
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APPENDIX A 

CLASY SYSTEM FLOWCHART 
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CLASy- 


-|SETUP9^ _ . . ^ 

rNXTCHR 

Lnumber 


•READTP 

rRINIT " 
■TAPHDR 



-FLDINT 

r-UNPACK 


.LINERD 

•UNPAKl 


-FOUNT 

-RREAD" 

-RWRITE- 

-ZOR 

L-BYTRAN+ 


.LAREAD- 

NXTCHR 


kMERR 


Imulti 1 

DATFIX 


* 1 

-ALFREE 

free;*; 

; i 

■CLINIT 

GET 

!' I 

’ 1 

-STATIS j 

1 DISC 


RREAD" 

CORRECT 

DOTSQ 

VPV 

VMTV 

MPVS 

ADJUST* 


STOFLO 

NRAND 


-CLDUMP rCLPR* 

Lisplit 


l-CLUSMP 


^BFINIT® 

^FDLINT 



rlSPLIT 
-CORRECT 
L DOTSQ 


+ = Univac Assembly Language Routine 
0 - Fortran internal subroutine 

= Univac Random I/O "random file access) routine 
* = Remainder of flowchart on page A- 3 
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ADJUST- 


-get— 

-TR 

-DOTSQ 

-AMSQ 

^SQMTX 

r-MINV 

5-UNIF- 

-CLPR 

-TRIMTX 

S-DENCAL 

-SPLIT- 


-FREE 

-CLDUMP- 

‘“SEPER-- 


rSUBLIM- 
-ELIM - 

-CORECT 

•JOIN 


■STOFLO 


•NRAND 


SQMTX 
EIGROT' 
ACOM - 
fMLT 
MVEC 
hMTVEC 
ACOM 
pOET- — 
l-APRIOR 
MINV 

CLPR 


I 


ISPLIT 
^CLPR 
1- DEN CAL 
^FREE 
-TRFREE. 
-rSUBLIM- 

Ltrfree- 

GET 

SQMTX 
MINV 
APRIOR 
DENCAL 
KRIMTX 
^CLPR 


FREE 
GET - 


■STOFLO 


STOFLO 


■ FREE 
■TRFREE. 
FREE 

STOFLO 


■ TREE 


l-APRIOR 




CLPR 


l-LOCK 
-GET- 
-SQMTX 
-MINV 
•-FREE 


STOFLO 
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